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ABSTRACT 

We explore the effects on nucleosynthesis in Type II supernovae of various parameters (mass cut, 
neutron excess, explosion energy, progenitor mass) in order to explain the observed trends of the iron- 
peak element abundance ratios ([Cr/Fe], [Mn/Fe], [Co/Fe] and [Ni/Fe]) in halo stars as a function of 
metallicity for the range —4 < [Fe/H] < —2.5. [Cr/Fe] and [Mn/Fe] decrease with decreasing [Fe/H], 
while [Co/Fe] behaves the opposite way and increases. We show that such a behavior can be explained 
by a variation of mass cuts in Type II supernovae as a function of progenitor mass, which provides a 
changing mix of nucleosynthesis from an alpha-rich freeze-out of Si-burning and incomplete Si-burning. 
This explanation is consistent with the amount of ejected 56 Ni determined from modeling the early light 
curves of individual supernovae. We also suggest that the ratio [H/Fe] of halo stars is mainly determined 
by the mass of interstellar hydrogen mixed with the ejecta of a single supernova which is larger for larger 
explosion energy and the larger Stromgren radius of the progenitor. 

Subject headings: Galaxy: evolution — Galaxy: halo — nucleosynthesis — stars: abundances — 
supernovae: general 



1. INTRODUCTION 

Very metal poor stars provide important clues to in- 
vestigate early Galactic chemical and dynamical evolu- 
tion because they contain valuable information about this 
time. We can learn how the Galaxy evolved chemi- 
cally (as well as dynamically) in its early phase from the 
observation of these stars. Therefore, a large number 
of observations and abundance analyses have been per- 
formed. Until recently, no significant abundance change 
with metallicity had been observed. However, recent high- 
resolution abundance surveys have discovered interesting 
trends of [Cr/Fe], [Mn/Fe] and [Co/Fe] with respect to 
[Fe/H] ([A/B] = log lo (A/B)-log lo (A/B) ). Both [Cr/Fe] 
and [Mn/Fe] decrease with decreasing metallicity from 
[Fe/H] = -2.4 to -4.0, while [Co/Fe] increases (Figure 
|) (McWilliam et al. 1995a, 1995b; McWilliam 1997; see 
also Ryan et al. 1996 and reference therein). 

Over the last few years, these trends have been the sub- 
ject of controversy. McWilliam et al. (1995b) discussed 
the time-delay mechanism, which produces the abundance 
trends by the time-delay between SNe II with different 
masses or between SNe I and SNe II. They argued that this 
mechanism works for several species such as a-elements 
but has difficulties with the above trends. McWilliam et al. 
(1995b) also pointed out that the Galactic halo formation 
model proposed by Searle & Zinn (1978) is at odds with the 
time-delay mechanism. Instead, they suggested that the 
metallicity-dependent supernova progenitor mass function 



could explain these trends, while still being consistent with 
Searle & Zinn (1978). Audouze & Silk (1995) interpreted 
the observations as evidence that the Galaxy was poorly 
mixed and proposed that at most two or three supernovae 
could have contaminated any particular cloud. In favor of 
this hypothesis, Ryan et al. (1996) proposed that a varia- 
tion of the supernova explosion energy, which is connected 
with both the yields and the amount of ejecta dilution, 
could cause these trends. McWilliam (1997) and Searle 
& McWilliam (1998) suggested that diluting a primordial 
composition with solar composition supernova ejecta could 
explain the trends. In principle, hypothetical pre-galactic 
population III objects, such as very massive stars or pair 
creation supernovae (e.g., Bond et al. 1984), could be the 
cause of the observed behavior. However, none of these 
objects has been investigated in full detail in quantitative 
nucleosynthesis yield studies nor ever been incorporated 
in chemical evolution calculations. 

Timmes et al. (1995) constructed a Galactic evolution 
model using a grid of Type II supernova models (Woosley 
& Weaver 1995, hereafter WW95), but did not explain 
the above-mentioned abundance trends. In Timmes et 
al. (1995), only the region with [Fe/H] > —3.0 was in- 
vestigated, so that only limited information is available. 
However, their results seem to suggest that the metallicity 
effect alone cannot explain the observed trends. 

The purpose of this paper is to provide a qualitative 
explanation of the halo abundance data within the frame- 
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work of Type II supernovae (SNe II) explosive nucleosyn- 
thesis. It should be pointed out that self-consistent ex- 
plosion and nucleosynthesis calculations, which start from 
core collapse, have not yet successfully determined the 
mass cut between the central neutron star (or black hole) 
and the ejected envelope. Therefore, the mass cut of SNe 
II as a function of the progenitor mass is still an open 
question and can be treated as a free parameter. This in- 
vestigation focuses on the mass cut, which determines the 
amount of 56 Ni in SN II ejecta, and its possible connection 
with the amount of (stable) Cr, Mn, Co, and Ni ejected. 
The connection between stellar evolution, explosive yields, 
and galactic evolution, observed via abundances in low 
metallicity stars, is made by performing chemical evolu- 
tion studies in the framework previously described by Tsu- 
jimoto et al. (1995, 1997). 

Although we are investigating the nucleosynthesis 
of low-metallicity stars, we use solar metallicity pre- 
supernova models in this paper for the following reasons. 
Previously, only WW95 calculated nucleosynthesis of SNe 
II using progenitor models with low metallicities; however, 
their progenitor models have not yet been published. Their 
results suggest that although the yields of iron peak ele- 
ments depend on the progenitor metallicity, their variation 
in the range Z = (1CT 2 - KT 4 )Z is only about 20% or 
less. Since we are interested only in Z in this range, we 
neglect metallicity effects on progenitors and explore the 
dependence on other parameters, such as the stellar mass, 
mass cut, explosion energy, and neutron excess (= 1 — 2Y e , 
where Y e denotes the electron mole number). With this 
procedure, we can investigate the trends of the yields, i.e., 
how much the yields increase or decrease with these pa- 
rameters. 

In §||, we review the nucleosynthesis yields of SNe II (in- 
cluding their dependence on the mass cut, explosion en- 
ergy, and Y e of the innermost ejected matter). The results 
are applied in chemical evolution calculations in §^ in or- 
der to explore whether the observed abundance patterns in 
the metal-poor stars can be explained with our approach. 
Section ^ contains a more in-depth discussion, listing both 
the successes and the remaining open questions. 

2. NUCLEOSYNTHESIS IN TYPE II SUPERNOVAE 
2.1. Models 

Stars initially more massive than 8M & explode as SNe 
II at the end of their life, leaving neutron stars or black 
holes. When they explode, nucleosynthesis proceeds ex- 
plosively because of high temperatures and densities in 
the deep stellar interior. In order to investigate the explo- 
sive nucleosynthesis of SNe II, we solve nuclear reaction 
networks, together with hydrodynamical equations, and 
calculate the energy generation by explosive burning. 

Our calculations are performed in two steps. The first 
step is a hydrodynamical simulation of the SN II explosions 
with a small nuclear reaction network which contains only 
13 alpha nuclei ( 4 He, 12 C, 16 0, 20 Ne, 24 Mg, 28 Si, 32 S, 36 Ar, 
40 Ca, 44 Ti, 48 Cr, 52 Fe, and 56 Ni). The hydrodynamical 
simulations are carried out with a one dimensional PPM 
(piecewise parabolic method) code (Colella and Woodward 
1984). We generate a shock by depositing thermal energy 
below the mass cut that divides the ejecta and the collaps- 
ing core, and perform the induced nucleosynthesis calcula- 



tions following the procedure of our previous investigations 
(e.g. Hashimoto et al. 1989; Thielemann, Hashimoto, & 
Nomoto 1990; Thielemann, Nomoto, & Hashimoto 1996, 
hereafter TNH96). In the second step, at each mesh point 
of the hydrodynamical model, post-processing calculations 
are performed with an extended reaction network (Hix & 
Thielemann 1996), which contains 211 isotopes, and pro- 
vides precise total yields (also for minor abundances). 

In our calculation, the progenitor models are taken from 
Nomoto & Hashimoto (1988, hereafter NH88). Figure | 
shows the Y e profiles of the models (25M Q , 20M© and 
IZMq) in NH88. In this paper, Y e in the deep stel- 
lar interior is modified to take a constant value of y e dee P 
from M r = M cut to M r = 1.70M© (for the 25M Q star), 
M r = 1.64M Q (2OM star) and M r = 1.49M Q (13M© 
star). NH88's Y e distribution of the progenitor is adopted 
for the outer region above the boundaries shown in Figure 
|2[ Some of the adopted mass cuts are smaller than the Fe 
core masses obtained in the stellar evolution calculations 
(NH88). Note, however, that the Fe core mass, the Y e 
profiles, and Y e dee P are subject to uncertainties involved in 
the stellar evolution calculations, such as the treatment of 
convection, reactions rates (especially 12 C(a, 7) 16 0)), etc. 
Using the modified models, we analyze the dependences 
on the mass cut, Y e , etc, and obtain some constraints on 
these parameters. We also treat mass cuts and explo- 
sion energies independently, though these are physically 
related, because of uncertainties in the explosion mecha- 
nism. When analyzing the yields, we do not include grid 
zones inside the mass cut in the ejecta by assuming that 
materials in these regions fall back onto a neutron star or 
a black hole. The change in the final total energy caused 
by this treatment is negligible. 

In the following subsections, we show the dependence 
of SNe II yields on various parameters (mass cut, neutron 
excess, explosion energy, progenitor mass). Among them, 
the mass cut seems to be the most important parameter. 

2.2. Production of iron-group elements in Type II 
supernovae 

Figure shows the isotopic composition of the ejecta 
of the explosion of the 20M Q star as a typical SN II. 
Here three distinctive regions are seen. One is the in- 
nermost high temperature region. The outer boundary 
of this region is determined by the condition that maxi- 
mum temperatures of T > 5 x 10 9 K is attained behind 
the shock. In this region, nuclear statistical equilibrium 
(NSE) is achieved except for the slow triple-alpha process 
(a- rich freezeout). Thus, the nucleosynthesis depends only 
on Y e and the entropy, being dominated by iron group el- 
ements (Woosley, Arnett, & Clayton 1973; Thielemann, 
Hashimoto, & Nomoto 1990). The most abundant ele- 
ment in this region of complete Si-burning with alpha-rich 
freeze-out is 56 Ni, provided Y e > 0.49. The second region, 
which does not experience such high temperatures, under- 
goes incomplete Si-burning with only partial Si exhaus- 
tion; here the most abundant nucleus changes from 56 Ni 
to 28 Si. In the third region the temperatures are too low 
to produce any iron group elements, and only explosive 
O, Ne, or C-burning take place. Farther out in radius, 
regions are encountered where explosive nucleosynthesis 
barely affects the pre-explosive composition and matter as 
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produced during the quasi-static evolution is ejected es- 
sentially unchanged. 

In order to investigate the ratios [Cr/Fe], [Mn/Fe] and 
[Co/Fe] we look into the regions where these elements 
are produced. First of all, 56 Ni, which decays into the 
most abundant Fe isotope 56 Fe, is produced not only in 
the complete Si-burning region but also in the incomplete 
Si-burning region. 52 Fe and 55 Co, which decay into 52 Cr 
and 55 Mn, respectively, are mostly synthesized in the in- 
complete Si-burning region. 52 Fe is also synthesized in the 
complete Si-burning region, but not a large amount. (Note 
that the ordinate of Figure || is log-scaled.) On the con- 
trary, 59 Cu, which decays into 59 Co, is produced in the 
complete Si-burning region. Note that 55 Mn and 59 Co are 
the only stable isotopes of these elements and therefore, 
their abundances are identical with the Mn and Co ele- 
ment abundances. 52 Cr dominates the Cr element abun- 
dance by 84% (in solar composition). Thus, it is sufficient 
to take these isotopes into consideration when discussing 
the abundances of [Cr, Mn, Co/Fe]. 



2.3. Dependence on mass cuts 

The above discussion suggests that the choice of the 
mass cut can affect the ratios [Cr/Fe], [Mn/Fe], and 
[Co/Fe]. For a deeper mass cut (i.e. smaller M cut ), the 
ejected mass of the complete Si-burning region is larger 
(i.e., the masses of Fe and Co are larger), while the ejected 
mass of the incomplete Si-burning region remains the same 
(i.e., the masses of Cr, Mn, and Fe are the same), accord- 
ingly the ratios of [Cr/Fe] and [Mn/Fe] are smaller and 
[Co/Fe] is larger. For a mass cut at larger radii (larger 
M cu t) these ratios show the opposite tendency. Therefore, 
specific choices of mass cuts in SNe II might explain the 
behavior of [Cr/Fe], [Mn/Fe], and [Co/Fe] in the metallic- 
ity range -4 < [Fe/H] < -2.5. 

Table |l| summarizes the yields of the 25, 20 M© stars for 
various mass cuts. Figure]! shows the dependence of the 
abundance ratios on the mass cut. Here we use the solar 
abundances by Anders and Grevesse (1989). For smaller 
■^cut, [Cr/Fe] and [Mn/Fe] are smaller, while [Co/Fe] and 
[Ni/Fe] are larger. Note that stable Ni is dominated by 
fe8 Ni and, thus, [Ni/Fe] is dominated by the 58 Ni/ 56 Ni ra- 
tio. One finds substantial 58 Ni production only in the 
inner complete Si-burning region, so that the behavior of 
[Ni/Fe] is similar to [Co/Fe]. The observational trends in 
Cr, Mn, Co and Ni are reproduced simultaneously. 



2.4. Dependence on neutron excess 

Figure || shows the Y e dependence of SN II yields for 
ydeep = q.4985 (left) and 0.4950 (right). The progenitor 
is the 25M star with the 8M He core and the explosion 
energy is 1.0 x 10 51 ergs. Table || summarizes the nucle- 
osynthesis products and compares the abundance ratios 
with the solar ratios for different values of y e deep . Figure^ 
shows how these ratios depend on Y e . For smaller Y e (i.e., 
more neutron-rich environment), more 59 Cu (which decays 
into 59 Co) and 58 Ni are produced, while the yield of 52 Fe 
is smaller. 58 Ni is more sensitive to Y e than 59 Co. The 
56 Ni yield is smaller for smaller Y e . Thus, all the metal to 
iron ratios are larger for smaller Y e . 



2.5. Dependence on explosion energy 

The explosion energy of a typical SN II is considered to 
be E ~ 1 x 10 51 ergs. For example, the explosion energy of 
SN 1987A has been estimated to be E = (1.0 - 1.5) x 10 51 
ergs from modeling of the early light curve (Shigeyama et 
al. 1987, 1988, 1990; Woosley et al. 1988; Nomoto et al. 
1997; Nakamura et al. 1998). Figure ^ shows the energy 
dependence of Type II supernova yields for the same pro- 
genitor model, i.e., the 8M He core of the 25M star. 
Here Y" e docp = 0.4985, while E = 0.4 x 10 51 ergs (left) and 
2.0 x 10 51 ergs (right). 

It is seen that the larger explosion energy forms a larger 
region of complete silicon burning, while the incomplete 
silicon burning region is enlarged only slightly. In other 
words, for larger E, the region of incomplete silicon burn- 
ing is shifted outward in mass without an appreciable 
change of the enclosed mass. Therefore, the SN II with 
a larger explosion energy produces a larger amount of iron 
group elements if mass cuts are the same. That is, a larger 
explosion energy affects abundance ratios in the same way 
as a smaller mass cut. The explosion energy of SNe II 
can vary by a factor of five (Burrows 1998). We find that 
changing explosion energies from 0.4 x 10 51 to 2.0 x 10 51 
ergs causes some changes in the abundance ratios, but can- 
not explain the large variations in the observations. 

2.6. Dependence on progenitor mass 

Figure |s| shows the progenitor mass dependence of SN 
II nucleosynthesis. The progenitors are the 6M He core 
of the 2OM star (left) and the 8M He core of the 25M 
star (right). Both models have the same kinetic energy of 
the ejecta (E = 1 x 10 51 ergs), but the deposited ener- 
gies are different because of the difference in the gravita- 
tional potential energies between the two models. Y^ eep 
is set to be 0.4985. Note that a more massive supernova 
has more massive complete and incomplete Si burning re- 
gions. As one can see in Figure ||, the mass ratio between 
the complete and incomplete Si-burning layers depends on 
the progenitor mass and the mass cut. 

3. ABUNDANCE RATIOS IN METAL POOR STARS 

More massive stars evolve faster. Thus, we expect the 
ejecta of the most massive stars to dominate the earliest 
phase of Galactic evolution, i.e. the period correspond- 
ing to the lowest [Fe/H]. If the mass cut M cut between 
the ejecta and the neutron star tends to be smaller for 
the larger mass progenitor which provides earlier chem- 
ical input into the galactic evolution, one could expect 
to reproduce the observed trend in [Cr/Fe], [Mn/Fe] and 
[Co/Fe]. Whether this trend agrees quantitatively with 
the observations has to be tested. 

The above trend between the progenitor mass and the 
mass cut cannot be monotonic because of two effects which 
compete with each other. One is the amount of neutrino 
absorbing matter, and the other is the depth of the gravi- 
tational potential. For intermediate massive stars, heavier 
stars eject relatively larger amounts of 56 Ni because of 
larger neutrino absorbing regions, while in more massive 
stars, such as SN1997D, the deeper gravitational poten- 
tial wins and 56 Ni is scarcely ejected due to the fall back 
(Burrows 1998). Note that recent observations, relating 
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the total 56 Ni mass to progenitor mass, seem to support 
this idea (see Figure |l4|). 

We expect that the trend of increasing 56 Ni yield with 
increasing progenitor mass holds over a large range of 
progenitor masses, while the reverse trend occurs more 
abruptly. Therefore, we focus here on the increasing trend 
(decreasing mass cut with increasing progenitor mass) up 
to a maximum progenitor mass, above which we assume 
that SNc II do not contribute to galactic chemical evolu- 
tion. 

Figure ^| shows the abundance ratios of various SN II 
models given in Table ^, where their parameters (mass 
cut, etc) are summarized. The abscissa of Figure ^ in- 
dicates only a model sequence (A-I). For these models, 
we take into account the observational trend between the 
ejected 5f> Ni mass and stellar mass by assuming that the 
5(> Ni yield of the 25Mq star is larger than that of less mas- 
sive stars (see Figure [u]). Correspondingly, Figure ^| shows 
that deeper mass cuts yield smaller [Cr/Fe] and [Mn/Fe], 
while they lead to larger [Co/Fe]. This means that if su- 
pernovae at smaller [Fe/H] have deeper mass cuts, the ob- 
served trends can be explained. 

We consider two models to relate [Fe/H] and the pro- 
genitor masses. One is the 'well-mixed' model, and the 
other is the 'unmixed' model. In the 'well-mixed' model, 
we assume that the mixing in the galactic scale is so effi- 
cient that the chemical uniformity is achieved in the early 
Galaxy. In the 'unmixed' model, on the contrary, we as- 
sume that the mixing is not effective and the composition 
of a metal-poor star is the same as produced by a single 
SN II. 

For these models, [Fe/H] is related to the stellar masses 
as follows. For the 'well-mixed' model, [Fe/H] at galactic 
time to, as well as other abundances (e.g. Cu, Mn, and 
Co), can be obtained by integrating the yields from all 
SNe exploded for t < to- [Fe/H] decreases monotonically 
towards the past and the age at the lower [Fe/H] corre- 
sponds to the shorter lifetime of more massive stars. In 
this way, we construct the Galactic halo model where the 
ages at [Fe/H] = —4, —3 and —2.5 roughly coincide with 
the time when stars of 25M Q , 20M Q and 13M Q exploded, 
respectively. Note that the mass vs. [Fe/H] relation is 
based on several assumptions (star formation rate, initial 
mass function, etc.), thus being subject to some uncer- 
tainties. The solid curves in Figure show the changes 
in the iron-peak abundance ratios as a function of [Fe/H]. 
To calculate these ratios, a set of models [B, G, H] is used 
in Figure |, and the Salpeter IMF of <f>(M) cx M~ 2 - 35 is 
adopted. The observed trends of [Cr/Fe] and [Mn/Fe] are 
relatively well reproduced in Figure [h]. 

In the 'unmixed' model, the next generation stars were 
born in the ejecta of a SN II which is uncontaminated 
by previous supernovae. Then these stars have the same 
heavy element abundances as a single SN II. In this case, 
[Fe/H] depends on how much interstellar material is mixed 
with the ejecta before forming the next generation stars. 
Thus it is not trivial to see why there is a re lation between 
[Fe/H] and the progenitor mass M, In §4.1 we will discuss 
that the relation can be determined by Af-dependence of 
the explosion energy and the Stromgren radius. For now, 
we assign M = 25, 20, and 13 M to [Fe/H] = -4, -3, 
and —2.5, respectively, though the exact correspondence 



depends on the M -dependent explosion energy and the 
Stromgren radius. The filled circles in Figure |l(J show the 
'unmixed' case for the model set [B, G, H]. The trends 
in [Cr/Fe] and [Mn/Fe] are well reproduced with larger 
contrasts than in the 'well-mixed' model. 

For other set of models, the contrast in the abundance 
ratios can be larger than the above case. The solid lines 
and the filled circles in Figure pT| use the set [A, G, I], and 
the contrasts in [Cr/Fe] and [Mn/Fe] are larger than the 
set [B, G, H] for both the well-mixed and unmixed cases. 

Though it is hard to argue which model, between 'well- 
mixed' and 'unmixed', is in better agreement with the ob- 
servations, the slope of the abundance ratios for the 'well- 
mixed' model seems to be slightly too gentle, while the 
'unmixed' model can produce a steeper slope. Especially, 
the 'unmixed' model better reproduces the inclination of 
the [Co/Fe] vs. [Fe/Hl curve. 

In Figures |l0| and [ll], the calculated [Co/Fe] ratios are 
smaller than the observed ratios by a factor of 3-5. In 
order to compare the slope of the curve with observations 
more easily, the dashed lines in Figures |l^ and [ll] show the 
models in which the amount of Co is five times larger than 
in the models [B, G, H] and [A, G, I]. The filled triangles 
in Figures |l^ and [ll] indicate the five time s en hanced Co 
for the 'unmixed' model. As discussed in §2.4, [Co/Fe] is 
larger for smaller Y e . Among the models in Figure ^| and 
Table § 'C', 'D', and 'F' have the same mass cut as 'A', 
'B', and 'G', respectively, but have smaller Y e and thus 
larger Co. The solid lines for the 'well-mixed' model and 
the filled squares for the 'unmixed' model in Figures R2 
and [l3] show [D, G, H] and [C, F, I], respectively. Ah 
though [Co/Fe] is larger, stable Ni ( 58 Ni, 60 Ni, and 62 Ni) 
is overproduced and [Ni/Fe] is too large to be acceptable. 
Therefore, the deficiency of 59 Co (i.e., the deficiency of 
59 Cu) is still an open question. 

4. DISCUSSION 

4.1. Contamination by Type II supernovae in the very 
early Galaxy 

In §|l we have shown that the trends of [Cr/Fe] and 
[Mn/Fe] can be well explained with both the 'well-mixed' 
model and the 'unmixed' model. However, the 'well- 
mixed' model tends to predict too small a metallicity de- 
pendence in the abundances ratios. Indeed, Audouze & 
Silk (1995) and Ryan et al. (1996) argue that only one, 
or at most two to three, supernovae could have contami- 
nated a particular cloud in the [Fe/H] < —2.5 region. For 
the 'unmixed' model, the observed variations of [Cr/Fe], 
[Mn/Fe], and [Co/Fe] can be explained with the system- 
atic variations of these ratios as a function of the progen- 
itor mass M (§j3|). The question is then why there should 
be a tight correlation between [Fe/H] and M. 

The [Fe/H] of a star is determined by the amount of iron 
ejected from the relevant supernova and the mass of hy- 
drogen in the mixing region. Our calculations show that in 
order to explain the large variations of the observed abun- 
dance ratios (e.g. —0.5 dex in [Cr/Fe]), the iron mass from 
SN II varies within a relatively narrow range of 0.05M Q to 
0.25M©. Thus in the very low metallicity regions, there 
must be an order of magnitude variation in the hydrogen 
mass to produce the variation of [Fe/H] in the range of —4 
to -2.5. 
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Ryan et al. (1996) obtained an analytic expression for 
the mass of the ISM (interstellar matter), mjgM, mixed 
with the ejecta using equations (4.4b) and (3.33a) in Cioffi 
et al. (1988). Ryan et al. (1996) suggested that ttiism 
depends only weakly on the environmental details, but 
strongly on the explosion energy E of the supernova as 
miSM oc E 95 . In this connection, we should note the 
recent discovery of a supernova with such a large explo- 
sion energy as E ~ 3 x 10 52 erg from a massive progen- 
itor (SN1998bw; Iwamoto et al. 1998). If more massive 
supernovae tend to produce larger explosion energy and 
this M-dependence is large enough, the larger M - smaller 
[Fe/H] relation can be obtained. 

Here we propose another possibility. The above dis- 
cussion is derived under the assumption that the ISM is 
uniform. However, there is a distinctive non-uniformity 
characterized by the "Stromgren sphere" , within which 
matter is ionized. The progenitors of SNe II are so hot 
and luminous during their main-sequence phase that their 
Stromgren spheres are as large as 10 ~ lOOpc. The 
Stromgren sphere is likely to determine the amount of hy- 
drogen into which the ejecta of a supernova are mixed. The 
shock advances easily within the Stromgren sphere, but is 
strongly decelerated outside the sphere because it has to 
ionize the matter there; this means the effective adiabatic 
index 7 approaches 7 ~ 1 outside the sphere. The radius of 
the Stromgren sphere is sensitive to the effective temper- 
ature of the star. The change in the effective temperature 
by 25% leads to a doubling of the Stromgren radius (Os- 
terbrock, 1989), which results in a factor of 2 3 = 8 change 
in the swept up hydrogen mass. Therefore, it is possible 
to explain the observed range in [Fe/H]. Also the larger M 
- smaller [Fe/H] relation can be explained, because a more 
massive, hotter star has a larger Stromgren radius. 

This possibility can be corroborated by time scale esti- 
mates. The recombination time scale, i.e., the Stromgren 
sphere's life time, is 10 5 - 10 6 years. The stellar evolution 
time after the main-sequence phase, plus shock propaga- 
tion time, is 10 5 - 10 6 years. Therefore the Stromgren 
sphere produced during the main-sequence phase survives 
until the SN II explosion. Then the shock radius coin- 
cides the Stromgren radius, which depends on both the 
progenitor mass and the number density of hydrogen in 
the ambient ISM per cubic centimeter, no. We esti- 
mate that the hydrogen mass swept up by the shock is 
4 x 10 4 ■ (no/l.O) _1 M for M = 25M Q and 4 x 10 3 • 
(no/l.O) _1 M for M = 15M , using low metal progen- 
itor models (Schaller et al. 1992). These estimates show 
that the hydrogen masses differ by a factor of 10. How- 
ever, we should note that, for very small no (no ~ 0.1), 
the shock cannot reach the Stromgren radius, so that the 
hydrogen mass does not depend on the Stromgren radius. 
On the other hand, for larger no (no ~ 1), the Stromgren 
sphere vanishes before being approached by the shock. In 
this case, however, the difference in the swept up mass 
between larger and smaller mass stars is larger, because 
the longer life time of the latter makes it easier for the 
Stromgren sphere to vanish. This possibility is worth in- 
vestigating. 



4.2. Mass of 56 Ni in Type II supernova ejecta 



Our results in §|3j suggest that a supernova from more 
massive progenitor has a deeper mass cut, thus ejecting 
a larger amount of 56 Ni mass. The mass of 56 Ni should 
be determined by the competition between the amount 
of neutrino absorbing matter and the depth of the gravi- 
tational potential. Accordingly, the intermediate massive 
stars eject a relatively large amount of 56 Ni because of a 
large neutrino absorbing region, while in a more massive 
star the deeper gravitational potential wins and 56 Ni is 
scarcely ejected due to fallback. 

This is consistent with the recent 56 Ni mass estimates 
from the light curve modeling for SNe II and Type Ib/Ic 
supernovae (SNe Ib/Ic) which are the explosions of bare 
cores of massive stars (Figure |lj). From SN lib 1993 J 
(Nomoto et al. 1993; Iwamoto et al. 1997), SN Ic 19941 
(Nomoto et al. 1994a; Iwamoto et al. 1994) and SN1987A 
(e.g. Nomoto et al. 1994b for a recent review), the 13 - 
20 Mq stars eject ~ 0.07M Q 56 Ni. For M > 20M Q , there 
has been little information, but recent SNe II/Ic suggest 
the strong mass- dependence of Fe yield. The light curve 
of SN Ic 1997ef indicate the production of ~ O.15M of 
56 Ni (Iwamoto et al. 1998), while that of SN II 1997D 
shows the synthesis of only ~ 0.002M© of 56 Ni from the 
25 - 30M Q star (Turatto et al. 1998). SN II 1994W is 
also a small Fe producer (Sollerman et al. 1998). Such an 
observed mass dependence of the 56 Ni mass is consistent 
with the theoretical expectation. 

For M > 3OM , we expect that little 56 Ni is ejected 
(WW95), because of a large amount of fall back. On the 
other hand, the outer mantle including the oxygen-rich 
layer would be ejected from such a massive supernova, 
which is also required from the Galactic chemical evolution 
model (Tsujimoto et al. 1997). Possible contributions of 
recently discovered hypernovae (Iwamoto et al. 1998) are 
worth investigating. 

4.3. The Galactic chemical evolution 

In our new models, stars around ~ 25Mq have deeper 
mass cuts, thus ejecting a larger amount of Fe compared 
with the yields adopted by Tsujimoto et al. (1995). It 
would be interesting to examine the consequences of such 
a stellar-mass dependent Fe yield for the Galactic chemi- 
cal evolution model and to compare with observations. For 
[Fe/H] < —2.5, however, the error bars of the abundances 
in metal-poor stars other than Cr, Mn, Co, Ni and Fe 
are still too large to make meaningful comparisons (Ryan 
1998). 

For [Fe/H] > —2.5, we can assume that the ejecta from 
various SNe II are well mixed with ISM. Thus the abun- 
dance ratios between various elements and Fe are deter- 
mined by the total amount of Fe ejected by previous SNe 
II. With the new choice of mass cuts, the total amount of 
Fe is almost the same as that in the previous model by 
Tsujimoto et al. (1995), because Tsujimoto et al. (1995) 
assumed that 13-15 Mq stars produce a larger amount of 
Fe than stars more massive than 20Af©, while in our model 
13 - 2OM0 stars produce smaller amounts of Fe than more 
massive stars. These differences almost cancel out, leading 
to similar total Fe yield from SNe II. 

Because the total Fe yield from SNe II is similar to Tsu- 
jimoto et al. (1995), the evolution of the abundance ratios 
such as [Cr/Fe] and [O/Fe] at [Fe/H] > -2.5 in our model 



6 



Nucleosynthesis in Type II supernovae and the abundances in metal-poor stars 



is similar to those in Tsujimoto et al. (1995). Figure |l5] 
shows the evolution of [O/Fe] against [Fe/H] for the differ- 
ent Fe yields from SNe II, i.e., our model (solid line) and 
Tsujimoto et al. (1995; dashed line). In our model, we use 
the set [B, G, H] and assume no Fe from stars more mas- 
sive than 26M0. The upper mass and the lower mass of 
IMF are set to be 40M Q and 0.1M Q , respectively, and the 
parameters of SNe la are chosen to reproduce the abun- 
dance distribution in the solar neighborhood (Kobayashi 
et al. 1998). These parameters are only slightly different 
from Tsujimoto et al. (1995). Our model and Tsujimoto et 
al. (1995) reproduce the observations almost equally well. 
The evolution of the [Mn/Fe] ratio between [Fe/H] = -2.5 
and —1.0 depends more on the SNe la rate rather than the 
SNe II rate (Kobayashi 1998). At [Fe/H] < -1, [Mn/Fe] 

0.5 due to SNe II. With increasing [Fe/H], [Mn/Fe] 

approaches the solar value because of the increasing con- 
tribution of SNe la, which have [Mn/Fe] > according 
to the SN la model W7 (Nomoto et al. 1984). Detailed 
Galactic halo chemical evolution models will be given in a 
separate paper. 

The deficiency of 59 Co (i.e., 59 Cu) remains a prob- 
lem. [Co/Fe] does not reach zero at [Fe/H] = in the 
above Galactic chemical evolution model. The most im- 
portant parameter which affects the abundance of 59 Cu 
is y e deop (see §[E|). Lower Y" e doc P produces larger [Ni/Fe] 
and [Co/Fe], which results in a large deviation of [Ni/Fe] 
from the observations. Therefore it seems unlikely that Y e 
solves this problem. Could uncertainties in nuclear cross 
sections solve this problem? This may not be easy because 
NSE is almost achieved in the region that produces 59 Cu. 
Thus, cross sections would not have a large effect. There 
are other elements which decay into Co, but the amount 
of these elements is very small. The deficiency of Co has 
to be investigated further. 

In this work, we use progenitor models with solar metal- 
licity. Metallicity dependence may not affect our main 
conclusions that the mass cut can explain the interest- 
ing trends in the abundance ratios, because the produc- 
tion sites of the relevant elements and their relative zone- 
thickness may not be so sensitive to metallicity. However, 
some quantitative features, e.g., the absolute value of each 
yield, may be sensitive to metallicity. Two of us (H. U. 
& K. N.) are currently investigating such metallicity ef- 
fects. There are other potentially important effects, such 
as instabilities, dynamic convection, and rotation. These 
effects also should be investigated in the future. 
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Fig. 1. — The trends of iron-peak abundance ratios. These data are from Ryan et al. (1991, 1996) (filled circles); 
McWilliam et al. (1995a), McWilliam (1998) (filled pentagons); Norris et al. (1993), Primas et al. (1994) (filled triangles); 
Gratton & Sneden (1987, 1988, 1991), and Gratton (1989) (open circles); Magain (1989), and Zhao & Magain (1990) 
(open pentagons); and Molaro & Bonifacio (1990), Molaro & Castelli (1990), Peterson et al. (1990) (open triangles). 
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Fig. 2. — Y e profiles of the models in NH88. Dotted lines indicate the boundaries of Y e modification. In our calculations, 
Y e is modified to be a constant value below these lines, while the NH88's Y e distribution of the progenitor models is 
maintained above these lines. 
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Fig. 3.— The isotopic composition of the ejccta from 20M Q stars (6M© He core). 
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Fig. 4. — Dependence on mass cut. Three models with different mass cuts are used for each line. Solid line indicates 
supernovae of M coro = 6M Q , E cxp = 1.0 x 10 51 ergs, Y dccp = 0.4940. Dotted line indicates that of M corc = 8M , 
E cxp = 1.0 x 10 51 ergs, Y c dccp = 0.4985. 
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Fig. 8. — Dependence on progenitor mass. Both models have the same explosion energy of 1.0 x 10 51 ergs and a Y^ ccp of 
0.4985, but different progenitors, which are the 6M Q He core of the 20M Q star (left) and the 8M Q He core of the 25M Q 
star (right). 
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Fig. 9. — Abundance ratios of various SNe II. Mass cuts of the models are also shown below abundance ratios. Models 
A and C (B and D, F and G) have different Y" e dcop values (See Table §). 
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Fig. 10. — Variations in the iron-peak abundance ratios. The observational trend between ejected 56 Ni mass and stellar 
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Fig. 13. — Same as Figure [ll|, but with models for making stronger contrast in abundance ratios. The models [C, F, I] 
in Figure | and Table | are used. 
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Fig. 14. — Ejected 56 Ni mass as a function of main-sequence mass, as estimated for SN II 1987A (Shigeyama & Nomoto 
1990); SN lib 1993 J (Shigeyama et al. 1994); SN Ic 19941 (Iwamoto et al. 1994); SN II 1997D (Turatto et al. 1998); SN 
Ic 1997ef (Iwamoto et al. 1998). The 56 Ni mass of the models in Figure ^ are shown as open circles. 
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Fig. 15.— The evolution of [O/Fe] against [Fe/H] for the different Fe yields from SNe II, i.e., our model with the set 
[B, G, H] (solid line) and Tsujimoto et al. (1995; dashed line). The observational data are taken from Edvardsson et al. 
(1993) (filled circles); Barbuy & Erdelyi-Mendes (1989) (open triangles); Nissen et al. (1994) (stars); and Gratton (1991) 
(open squares). 
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Dependence on mass cut. 
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solar ratio 


Cr/Fc 


5.34E-03 


7.87E-03 


1.09E-02 


7.62E-03 


1.26E-02 


2.08E-02 


1.40E-02 


[Cr/Fe] 


-0.417 


-0.249 


-0.106 


-0.263 


-0.044 


0.173 




Mn/Fc 


1.38E-03 


2.08E-03 


2.96E-03 


1.95E-03 


3.50E-03 


6.05E-03 


1.04E-02 


[Mn/Fe] 


-0.878 


-0.700 


-0.548 


-0.728 


-0.475 


-0.273 




Co/Fe 


2.82E-03 


2.77E-03 


2.64E-03 


3.10E-03 


1.27E-03 


2.45E-03 


2.64E-03 


[Co/Fe] 


0.028 


0.022 


0.001 


0.070 


-0.032 


-0.316 




Ni/Fe 


3.57E-01 


3.39E-01 


3.20E-01 


3.74E-01 


2.83E-01 


1.30E-01 


5.76E-02 


[Ni/Fc] 


0.792 


0.770 


0.745 


0.812 


0.692 


0.355 





Table 2 
Dependences on Y e . 



Model 



M(M ) 


25 


25 


25 


20 


20 


20 




M corc (M Q ) 


8 


8 


8 


6 


6 


6 




E cxp (xl0 51 ergs) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 




M cut (M Q ) 


1.42 


1.42 


1.42 


1.42 


1.42 


1.42 




yrdccp 
* e 


0.4950 


0.4965 


0.4985 


0.4940 


0.4965 


0.4985 




Yield (M Q ) 














main clem. 


Fe 


2.44E-01 


2.71E-01 


3.05E-01 


1.72E-01 


1.99E-01 


2.21E-01 


b6 Ni 


Cr 


1.03E-03 


1.35E-03 


1.41E-03 


1.31E-03 


1.37E-03 


1.44E-03 


52 Fe 


Mn 


3.38E-04 


3.38E-04 


3.39E-04 


3.35E-04 


3.37E-04 


3.38E-04 


55 Co 


Co 


6.87E-04 


6.05E-04 


4.03E-04 


5.33E-04 


4.47E-04 


3.05E-04 


59 Cu 


Ni 


8.71E-02 


6.31E-02 


3.26E-02 


6.42E-02 


4.05E-02 


2.29E-02 


58 Ni 


Ratio 














solar ratio 


Cr/Fc 


5.34E-03 


4.97E-03 


4.61E-03 


7.62E-03 


6.89E-03 


6.51E-03 


1.40E-02 


[Cr/Fc] 


-0.417 


-0.448 


-0.481 


-0.263 


-0.306 


-0.331 




Mn/Fc 


1.38E-03 


1.25E-03 


1.11E-03 


1.95E-03 


1.69E-03 


1.53E-03 


1.04E-02 


[Mn/Fe] 


-0.878 


-0.922 


-0.973 


-0.728 


-0.790 


-0.833 




Co/Fe 


2.82E-03 


2.34E-03 


1.32E-03 


3.10E-03 


2.25E-03 


1.38E-03 


2.64E-03 


[Co/Fe] 


0.028 


-0.072 


-0.300 


0.070 


-0.070 


-0.281 




Ni/Fc 


3.57E-01 


2.33E-01 


1.07E-01 


3.74E-01 


2.03E-01 


1.04E-01 


5.76E-02 


[Ni/Fc] 


0.792 


0.607 


0.268 


0.812 


0.548 


0.256 
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Table 3 

Type II supernova Models in Figure ||. 



Model 


A 


B 


C 


D 


E 


F 


G 


B 


I 


M(M ) 


25 


25 


25 


25 


20 


20 


20 


13 


13 


M coro (M ) 


8 


8 


8 


8 


6 


6 


6 


3.3 


3.3 


E exp (xl0 51 ergs) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


M cut (M ) 


1.42 


1.47 


1.42 


1.47 


1.42 


1.58 


1.58 


1.32 


1.37 




0.4985 


0.4985 


0.4950 


0.4950 


0.4940 


0.4940 


0.4985 


0.4990 


0.4990 


Yield (M ) 


Fe 


3.05E-01 


2.62E-01 


2.44E-01 


2.09E-01 


1.72E-01 


7.77E-02 


9.01E-02 


6.87E-02 


3.50E-02 


Cr 


1.41E-03 


1.39E-03 


1.30E-03 


1.29E-03 


1.31E-03 


1.18E-03 


1.21E-03 


1.04E-03 


9.85E-04 


Mil 


3.40E-04 


3.39E-04 


3.38E-04 


3.37E-04 


3.35E-04 


3.33E-04 


3.34E-04 


3.15E-04 


3.14E-04 


Co 


4.03E-04 


3.43E-04 


6.85E-04 


5.89E-04 


5.33E-04 


1.63E-04 


1.03E-04 


4.97E-05 


1.09E-05 


Ni 


3.26E-02 


2.75E-02 


8.71E-02 


7.36E-02 


6.42E-02 


1.84E-02 


7.92E-03 


3.79E-03 


1.09E-03 


56 Ni 


2.90E-01 


2.48E-01 


2.25E-01 


1.93E-01 


1.60E-01 


7.04E-02 


8.35E-02 


6.41E-02 


3.17E-02 


K,atio 


[Cr/FeJ 


-0.481 


-0.420 


-0.418 


-0.355 


-0.263 


0.035 


-0.018 


0.033 


0.305 


[Mn/Fc] 


-0.973 


-0.907 


-0.878 


-0.812 


-0.728 


-0.387 


-0.450 


-0.357 


-0.065 


[Co/Fel 


-0.300 


-0.304 


0.027 


0.027 


0.070 


-0.100 


-0.361 


-0.562 


-0.927 


[Ni/Fo] 


0.268 


0.262 


0.792 


0.785 


0.812 


0.614 


0.184 


-0.019 


-0.267 



